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Abstract 

We study inclusive 6-jet and 66-dijet production at the CERN LHC invoking the hypothesis 
of gluon Reggeization in i-channel exchanges at high energy. The 6-jet cross section includes 
contributions from open 6-quark production and from 6-quark production via gluon-to-bottom- 
pair fragmentation. The transverse-momentum distributions of inclusive 6-jet production measured 
with the ATLAS detector at the CERN LHC in different rapidity ranges are calculated both within 
multi-Regge kinematics and quasi- multi-Regge kinematics. The 66-dijet cross-section is calculated 
within quasi- multi-Regge kinematics as a function of the dijet invariant mass M,-,-, the azimuthal 
angle between the two jets A(f> and the angular variable \. At the numerical calculation, we adopt 
the Kimber-Martin-Ryskin and Bliimlein prescriptions to derive unintegrated gluon distribution 
function of the proton from its collinear counterpart, for which we use the Martin-Roberts-Stirling- 
Thorne set. We find good agreement with measurements by the ATLAS and CMS Collaborations 
at the LHC at the hadronic cm. energy of VS = 7 TeV. 
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I. INTRODUCTION 



The study of fo-jet hadroproduction provides an important test of perturbative quantum 
chromodynamics (QCD) at high energies. The total collision energies, a/5* = 1.8 TeV and 
1.96 TeV in Tevatron runs I and II, respectively, and \/S = 7 TeV or 14 TeV at the LHC, 
sufficiently exceed the characteristic scale fi of the relevant hard processes, which is of order 
of 6-jet transverse momentum p?, i.e. we have Aqcd <C A* <C \/~S. In this high-energy regime, 
so called " Regge limit" , the contribution of partonic subprocesses involving t-channel parton 
(gluon or quark) exchanges to the production cross section can become dominant. Then the 
transverse momenta of the incoming partons and their off-shell properties can no longer be 
neglected, and we deal with "Reggeized" t-channel partons. These t— channel exchanges obey 
multi- Regge kinematics (MRK), when the particles produced in the collision are strongly 
separated in rapidity. If the same situation is realized with groups of particles, then quasi- 
multi-Regge kinematics (QMRK) is at work. In the case of 6-jet and 66-dijet inclusive 
production, this means the following: 6-jet (MRK) or 66-dijet (QMRK) is produced in the 
central region of rapidity, while other particles are produced with large modula of rapidities. 

The parton Reggeization approach [lj is based on the hypothesis of parton Reggeization 
in t— channel exchanges at high energy [2]. It was used for the description of a large number 
of hard processes at the modern hadron colliders and the obtained results confirm the as- 
sumption of a dominant role of MRK or QMRK production mechanisms at high energy. This 
approach was successfully applied to interpret the production of isolated jets [21], prompt 
photons ]4], diphotons ]5], charmed mesons heavy quarkonia (7 10] measured at the Fer- 
milab Tevatron, at the DESY HERA and at the CERN LHC. The theoretical background of 
a parton Reggeization approach is the effective quantum field theory implemented with the 
non-Abelian gauge-invariant action including fields of Reggeized gluons [2J and Reggeized 



quarks 



111 ], which was proposed by L. N. Lipatov in 1995 [12]. In this effective theory 



Reggeized partons interact with quarks and Yang-Mills gluons in a specific way. Recently, 



in Ref . 



131 ] . the Feynman rules for the effective theory of Reggeized gluons were derived for 



the induced and some important effective vertices. 

Usually it is suggested that MRK or QMRK production mechanism to be the dominant 
one only at small px values. Our recent study of isolated jet production at the Tevatron 
and LHC colliders, see Ref. [2j , demonstrated that the parton Reggeization approach can be 
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successfully used already in the range of xt = ^ < 0.1, or at the px Hz 300—400 GeV for the 
energy \J~S = 7 TeV at the LHC. This result motivates us to apply the parton Reggeization 
approach for the study of 6-jet and 66-dijet production in the kinematical range of transverse 
momentum 20 < pr < 400 GeV and rapidity \y\ < 2.1, as it was measured by the ATLAS 
Collaboration at the CERN LHC Q. 

The high-energy factorization scheme with the effective vertices for Reggeized gluons 
has been used earlier in Refs. H] for description of inclusive open 6-quark b- 



jet [18| and 66-dijet [19| production at the Tevatron collider. In this paper, we study in 
the same manner the inclusive 6-jet and 66-dijet production at the CERN LHC invoking 
the hypothesis of gluon Reggeization in t-channel exchanges at high energy. We take into 
account two mechanisms of 6-jet production: the open 6-quark production and "jet-like" 



fe-quark production via gluon-to-bottom-pair fragmentation 20|]. We consider 6-quark jet as 



an isolated, by the jet-cone condition |2l|, hadronic jet containing one 6(6)-quark or 66-quark 
pair. Thus, the 6-jet production cross section can be written as a sum of two terms. The 
first one represents a so-called "open 6-quark" production, when the 6-jet contains 6(fe)-quark 
which is produced directly in the hard partonic subprocess. The second term corresponds to 
the case of "jet-like" production, where a 6-jet contains 66-quark pair which is produced via 
gluon or light-quark fragmentation. The transverse-momentum distributions of inclusive 
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in the different 



6-jet production measured with the ATLAS detector at CERN LHC 
rapidity ranges are calculated both within multi-Regge kinematics and quasi-multi-Regge 
kinematics. The 66-dijet cross sections are calculated within quasi-multi-Regge kinematics 
as functions of the 66-dijet invariant mass Mjj, the azimuthal angle between the two jets 
A0 and the angular variable x- 

This paper is organized as follows. In Sec. |TT]the parton Reggeization approach is briefly 
reviewed. We write down the relevant for our analysis analytical formulas for squared matrix 
elements and differential cross sections. In Sec. IHIt we describe our calculations and present 
the results obtained. In Sec. IIVI the conclusions are summarized. 

II. MODEL 

We study 6-jet production in the region of large 6-quark transverse momentum 3> Tn^, 
where mj is a 6-quark mass. At the present time, the conventional approach for calculation of 
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the b— quark production cross sections is based on the next-to-leading (NLO) approximation 



in perturbative QCD and collinear parton model [22|]. It is well known that fixed-order 
perturbation QCD calculations are applicable when the transverse momentum px of the 
produced heavy h— quark is not much larger than its mass ra&. In the case when the transverse 
momentum significantly exceeds the mass, the large logarithms of type log(pr/m&) arise to 
all orders of a s (//), so that a fixed-order approach breaks down [23|]. It is possible to resum 
all these logarithms in the fragmentation approach using the factorization theorem, which 
states that the cross section for the production process of high-p^ b— quark can be written 
in factorized form as a convolution of the short- distance partonic cross section of parton / 
production with the fragmentation function D b f{z, /i 2 ) for a formation of a b— quark from the 
parton /: 

Sr ^jdzJdp' T -^- r D b f (z,fj l )S(p T -zp' T ). (1) 



dpT j J J a P T 

The fragmentation functions for heavy quarks in perturbative QCD have been studied at 
the next-to-leading order (NLO) QCD approach in Ref. {24]. 

The experimentally measured transverse energy Et (or the transverse momentum p?) of 
6-jet includes transverse energies (transverse momenta) of all partons inside some jet-cone 



in the rapidity-azimuthal angle plane, which radius is defined as follows, R = a/ Ay 2 + A0 2 
2l\ . Such a way, it is insignificant which part of the initial parton four-momentum is 
transferred to the 6-quark, and we can simplify the formula ([1]) to the form 

fofrag fof 

— A — = Er n /W' 2 

dpr ^ dp T 

where nf(fji) = J D b (z, fj,)dz is a 6-quark multiplicity in the /-parton jet. It is obvious that 
a 6-quark multiplicity in a gluon-initiated jet greatly exceeds a 6-quark multiplicity in any 
quark- initiated jets, n g (fi) 3> n q ([i>) with q = u,d, s, c. We will take into account only main 
contribution from the gluon-to-bottom-pair fragmentation g — > bb. Let us note that in this 
case the 66-pair is considered as a one 6-quark jet. 

To describe inclusive 6-jet and 66-jet cross sections in terms of the parton Reggeization 
approach, in the LO we need to consider gluon fusion subprocesses of open 6-quark and 
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gluon production only, which are to be dominant at the high energy, they write: 

R(qi) + R(q2) -> g(p), (3) 
R(qi) + R(q 2 ) 6(pi) + 6(P2), (4) 
R(qi) + R(q 2 ) -> </(Pi)+^(P2), (5) 

where i? is a Reggeized gluon and g is a Yang-Mills gluon, respectively, with four-momenta 
indicated in parentheses. The contribution of the partonic subprocess (j3J) can be neglected 
in comparison with the contribution of the subprocess (jlj) because of the strong suppression 
by the g — >■ bb fragmentation (n g ~ 1CT 3 ) for both produced gluons. In Ref. [l6| it was 
shown, that at the Tevatron energy range, the contribution of the subprocesses Q + Q — >■ g 
and Q + Q — > & + 6 with initial Reggeized quarks is sufficiently smaller comparing to the 
dominant contribution of the subprocesses ([3]) and (HI), and the former becomes sizeable only 
at the very large 6-jet transverse momentum p?. As the LHC energy exceeds by a factor 3.5 
the one of the Tevatron collider, we estimate a quark-antiquark annihilation contribution to 
be even much smaller and therefore do not consider it in the present analysis. 

Performing a study of high-transverse- momentum b— quark production (jp? » m^) in 
the collinear parton model, we have an additional b— quark production mechanism, namely 
a production via 6— flavor excitation, where 6(6)— quarks are considered as partons in the 
colliding protons. For example, this mechanism has been used successfully to describe 
B— meson p?— spectra at the Tevatron and LHC in NLO calculations of the parton model 
We have used a similar idea in our previous study of inclusive b— jet production at the 
Tevatron within the Parton Reggeization Approach [16J. In this work we took into account 
the LO in a s contribution from 2 — > 1 partonic subprocess 

BR -»■ b, (6) 



where B is a Reggeized b— quark. As it is shown in the Fig. 1 of Ref. [161 ]. the sum of 
this contribution and a contribution from the subprocess (jlj) strongly overestimates the 
experimental data. In the present analysis we ignore a contribution from the subprocess 
(EJ). First, we avoid any chance of a double-counting between subprocesses (JH]) and (EJ). 
Second, the conception of quark Reggeization for a b— quark inside a proton seems to be 
wrong, b— quarks are produced preferably at the last step of QCD-evolution at the large 
scale [A ~ Pt, and their PDF is proportional to a large logarithm log(2>r/ra&). However, the 



QCD-evolution of a Reggeized parton should be valence- like. It means, the Reggeized parton 
must be a t— channel parton throughout all steps of QCD-evolution in the parton ladder. 
But, a 6— quark conventional collinear PDF, which we take as input for a KMR (Bliimlein) 
prescription to obtain a b— quark unintegrated PDF, satisfies sea-like QCD-evolution. For 
this reason we strongly overestimate a value of a b— quark unintegrated PDF. The more 
adequate way should be to consider a subprocess bR — >■ b with a collinear b— quark in 
the initial state, instead of subprocess ([6]). But even in this case a problem of double- 
counting still exist. That is why in the present study we consider Reggeized-gluon induced 
contributions, like ([3]) and (TjJ, only. 

The squared amplitude of subprocess ([3]) reads 



\M(R+R^g)\ 2 = -7ra s p 2 T , 



where p\ = t% + t 2 + 2^t x t 2 cos0i 2 , h 



<f 1T , t 2 



-Q2 



(7) 

q 2 T, with qiT and q 2 T 



representing the transverse momenta of initial Reggeized gluons, and <pi 2 is the azimuthal 
angle enclosed between them. 



The squared amplitude of subprocess (TjJ was obtained in Ref. 



291 ] using the effective 



Feynman rules of the parton Reggeization approach. It coincides with previous result of 



Ref. [30| which is expressed in the alternative form. The answer of Refs. 



29 



30t ] can be 



written down as a linear combination of an Abelian and a non-Abelian term, as 



\M(R + R-+ b + b)\< 



256vr 2 a 2 



2N, 



M A + 



2(N? - 1 



-M 
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t = t — ml, u = u — ml, a.\ 
(3 2 = 2{p 2 - Pi) / S . Here, the Mandelstam variables are defined as s = (qi + q2) 2 , t = (qi—pi) 2 , 
u = (^2 — Pi) 2 , S = (Pi + P2) 2 , where P\ and P2 denote the four- momenta of the incoming 
protons. 



7?i + (3 2 " a x + a 2 
2(p! ■ P 2 )/S, a 2 = 2(p 2 ■ P 2 )/S, fa = 2(pi ■ P 1 )/S, and 



Exploiting the hypothesis of high-energy factorization, we express the hadronic cross 
sections da as convolutions of partonic cross sections da with unintegrated PDFs $^ of 
Reggeized gluon in the hadrons h. For the processes under consideration here, we have 

d 2 q 2T 



da( PP ^gX) = [ *E± f [ *?1 [ 

J Xi J 7T J X 2 J 
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x *i, /i 2 )^(x 2 , t 2 , n 2 )da{RR ->■ g), 



daipp^bbX) = [ 

J Xi J 7T J X 2 J 



d 2 q 2T 



x m Xl , ti, ^ 2 )$^(a; 2 , * 2 , ^ 2 )da{RR ->■ 66). 



(10) 



(11) 



The unintegrated PDFs t,/i 2 ) are related to their collinear counterparts F^(x,fi 2 ) by 



the normalization condition 



rft^(x,t,/i 2 ), 



(12) 



which yields the correct transition from formulas in the parton Reggeization approach to 
those in the collinear parton model, where the transverse momenta of the partons are ne- 
glected. 

In our numerical analysis, we adopt as our default the prescription proposed by Kimber, 
Martin, and Ryskin (KMR) 3lJ to obtain unintegrated gluon PDF of the proton from the 
conventional integrated one, as implemented in Watt's code [321. The precise analysis of 
KMR gluon unintegrated PDF had been performed in the Ref. 33j, including an accurate 
study of the dependence on the choice of collinear input. As is well known 34| . other popular 
prescriptions, such as those by Blumlein or by Jung and Salam 36 1, produce uninte- 
grated PDFs with distinctly different t dependences. In our analysis we don't evaluate the 
unintegrated gluon PDF after Jung and Salam |36[] because this PDF had been tabulated 
only in a range of t, n 2 < 10 4 GeV 2 . It is not enough to calculate b — jet production cross 
sections up to pt = 400 GeV, in accordance with measurements of the relevant experiments. 
In fact, we had to use the unintegrated gluon PDF up to t,/i 2 < 10 6 GeV 2 . In order to 
assess the resulting theoretical uncertainty, we also evaluate the unintegrated gluon PDF 
using the Blumlein approach, which resums small-x effects according to the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) equation [2]. As input for these procedures, we use the LO set of 
the Martin-Roberts-Stirling-Thorne (MRST) [37] proton PDF as our default. The rele- 
vant theoretical study of particle production in the high-energy factorization scheme using 
KMR and Blumlein unintegrated gluon PDFs j^l-Q demonstrates that both unintegrated 



PDFs lead to a similar behavior of production spectra at the non-large particle transverse 
momentum (jp? < 20 GeV). In case of high transverse momentum production of isolated jets 
and prompt photons [3|, the theoretical predictions obtained with these PDFs are different. 
Although we take identical collinear inputs for both KMR and Bliimlein approaches, the 
relevant kernels of integrand transformation between collinear and unintegrated PDFs differ. 
The KMR approach is based on DGLAP evolution equation, while the Bliimlein approach 
is based on BFKL evolution equation. As the BFKL approach seems to be preferable in 
the region of very small x <C 1, which corresponds to non-large Pt at fixed \fS, the KMR 
unintegrated gluon PDF should be more suitable to describe the experimental data at large 
Pt- 

Throughout our analysis the renormalization and factorization scales are identified and 
chosen to be \i — £pt, where £ is varied between 1/2 and 2 about its default value 1 to 
estimate the theoretical uncertainty due to the freedom in the choice of scales. The resulting 
errors are indicated as shaded bands in the figures. 

The master formula for the doubly differential cross section of inclusive b— jet production 
via gluon-to-bottom-pair fragmentation at the pt 3> rra& reads as follows 

da fra9 (pp -> bX) 1 



dpT<iy p^ 



xn g (n)\M(RR^g)\ 2 , (13) 

where y is the rapidity of b— quark, 0i is the azimuthal angle enclosed between the vectors 
qi T and p T , 

PTexp(±y) 2 n /— , , s 

x 1;2 = -7= , h = h +Pt ~ 2y/t 1 p T cos{(p 1 ). 

In case of bb— dijet production via the partonic subprocess (HJ we get the differential cross 
section in the form: 



da open (pp -»■ bbX) p lT p 2T f 7 f 1± 2NjRO . \M(RR^bb)\ 2 

dp\Tdy\dp 2 Tdy 2 dl\(p \bn 6 J J [xix 2 b) z 

where Pi >2 t and y\ )2 are b— quark and b— antiquark transverse momenta and rapidities, re- 
spectively, A0 is the azimuthal angle enclosed between the vectors pyt and p 2 T, 

xi = (p? +p° 2 +pI + ti)/Vs, x 2 = (rf + pi - v \ - P z 2 )/^s, 

Pi,2 = t ex P(l/i,2) + exp(-yi >2 )] , p*a = ^%xp(y 1)2 ) - exp(-yi, 2 )]. 
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The inclusive 6-jet transverse-momentum spectrum can be presented in the form: 

rj n b i et rlrrf ra 9 rlrr°P en drr°P en 

5— = %— + 2 x ^— 6>(% -R) + ^9(R - Ra), (15) 
apx apx dpT dpx 



where R b i = a/ (yt — Ui) 2 + {4>b — 0b) 2 , R is the experimentally fixed jet radius parameter, 
6(x) is the unit step function. In such a way, the subprocess (00) of open 6-quark production 
contributes two separate 6-quark jets while R b i > R, and only one 6-quark jet while R b i < R. 



III. RESULTS 



Recently, the ATLAS Collaboration presented data on inclusive and dijet production 
cross sections which have been measured for jets containing 6-hadrons (6-jets) in proton- 
proton collisions at a center-of-mass energy of \J~S = 7 TeV |l4|. The inclusive 6-jet cross 
section was measured as a function of transverse momentum in the range 20 < pr < 400 
GeV and rapidity in the range \y\ < 2.1. The 66-dijet cross section was measured as a 
function of the dijet invariant mass in the range 110 < Mjj < 760 GeV, the azimuthal angle 
difference between the two jets A0 and the angular variable x i n two dijet mass regions. 
Jets were reconstructed with jet radius parameter R = 0.4. The angular variable \ is 
defined as follows x — ex P \yi — 2/2 1 - To measure the cross section as a function of x, an 
additional acceptance requirement was used that restricts the boost of the dijet system to 

IVboostl = + 2/2 1 < 1-1- 

The 66-dijet cross-section as a function of dijet invariant mass Mjj for 6-jets with 
Pt > 40 GeV and \y\ < 2.1 is shown in Fig. [lj The data are compared to LO parton 
Reggeization approach predictions, the solid polyline corresponds to KMR unintegrated 

n 1 1 

PDF [31|, the dashed one — to Bliimlein PDF [35]. We observe nice agreement between data 
and theoretical prediction obtained with the KMR unintegrated PDF. In case of Bliimlein 
PDF, the theoretical histogram lies about factor 2 lower than the experimental data and 
this difference increases towards the high values of dijet invariant mass. 

In Fig. |2] the 66-dijet cross-section as a function of the azimuthal angle difference A<p 
between the two jets for 6-jets with pt > 40 GeV, \y\ < 2.1 and a dijet invariant mass of 
Mjj > 110 GeV is presented. The normalized to the total cross section data are compared 
to LO parton Reggeization approach predictions, the solid polyline corresponds to KMR 
unintegrated PDF, the dashed — to Bliimlein PDF. For the both unintegrated PDFs our 
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predictions lie within the experimental uncertainty interval of data except only one point at 



the A(f> ~ 2. We need to mention that in the case of CDF measurements at the Tevatron [19 1 
the azimuthal-separation-angle distribution of inclusive 66-dijet production is well described 
using the parton Reggeization approach formalism at the all values of the azimuthal angle 
difference < A0 < tt (see Fig. 4 in the Ref. 16j). 

The 66-dijet cross-section as a function of angular variable x f° r 6-jets with px > 40 GeV, 
\y\ < 2.1 and \yboost\ = 2I2/1+Z/2I < 1-1, for dijet invariant mass ranges 110 < Mjj < 370 GeV 
and 370 < Mjj < 850 GeV are shown in the Fig. [3] and Fig. HJ correspondingly. The 
normalized to the total cross section data are compared to our LO parton Reggeization 
approach predictions. In the range of 110 < Mjj < 370 GeV the polylines corresponding 
to KMR and to Bliimlein unintegrated PDFs coincide. In the region of large invariant 
masses 370 < Mjj < 850 GeV, the prediction obtained with the Bliimlein unintegrated PDF 
lies about factor 2 lower than the data. On the contrary, the calculations with the KMR 
unintegrated gluon PDF are found to be in a good agreement with data. 

To calculate inclusive 6-jet transverse-momentum production spectra we need to take 
into account gluon-to-bottom-pair production mechanism and to use the function of 66- 
pair multiplicity n g (n) in a gluon jet. Because the existing theoretical predictions (see, for 
example, Ref. 25j) contain large uncertainties, we consider n 9 (/i) as a free phenomenological 
parameter, which is extracted from the experimental data from the ATLAS Collaboration 



14| for the inclusive 6-jet cross sections. 

In the Fig. El the inclusive differential 6-jet cross section as a function of pt for 6-jet s 
with \y\ < 2.1 is compared with our LO predictions of the parton Reggeization approach. 
The contribution of QMRK subprocess P| and the contribution of MRK subprocess ([3]) are 
shown separately. We see that open 6-quark production mechanism does not describe data, 
especially at the large pt and some contribution from gluon-to-bottom-pair fragmentation 
mechanism is needed. We have obtained good description of the data using n g (fi) as a 
free parameter. In Fig. El the 66-pair multiplicity n g (fi) in a gluon jet as a function of pt 
extracted from the ATLAS data for the inclusive 6-jet production spectra 14J is shown. The 
open circles and dashed fitting line correspond to Bliimlein unintegrated PDF, the black 
circles and solid fitting line correspond to KMR unintegrated PDF. The general theoretical 



consideration 



25] leads to the following analytical approximation for the 66-pair multiplicity 
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in a gluon jet 



^2 

rig(fi) = Aha.—, (16) 



2 ' 



where we fixed m h = 4.75 GeV and \i = px/4= |26|, and found that A KMR = 0.0012 in case 
of KMR unintegrated PDF, and As = 0.0027 in case of Bliimlein unintegrated PDF. At 
the scale [i ^ mz/4, which corresponds gluon-to-bottom-pair fragmentation of secondary 
gluon in the Z-boson decay [Z — > qq — > qqbb), our approximation ( I16p yields n g ~ 0.002 — 
0.004, that is in an agreement with the measurements at the LEP Collider: n g = (3.3 ± 
1.8) x 10~ 3 from the DELPHI Collaboration [38], n g = (2.44 ± 0.93) x 10" 3 from the SLD 



Collaboration 



391 ] . The difference in obtained 66-pair multiplicities n g (fi) with the KMR 



and Bliimlein unintegrated PDFs should be used to distinguish last ones. We conclude 



that KMR unintegrated PDF looks preferab 
Opposite this conclusion, we found recently 



v to describe 6-jet production cross sections. 
3[ that Bliimlein unintegrated PDF is better 



to describe inclusive all-flavor inclusive jet production spectra 



28] 



The measured by ATLAS Collaboration 14( inclusive double- differential 6-jet cross- 
sections as functions of for the different rapidity ranges: (1) \y\ < 0.3 ( x 10 6 ), (2) 
0.3 < \y\ < 0.8 ( x 10 4 ), (3) 0.8 < \y\ < 1.2 ( x 10 2 ) and (4) 1.2 < \y\ < 2.1 are shown 
in Fig. [71 Here, our theoretical predictions are obtained taking into account both contri- 
butions, open 6-quark production and fragmentation production with n g (fi) as in (fl6l) . and 
with KMR unintegrated PDF. We demonstrate good agreement with data in all rapidity 
intervals. 

To test the universality of the approach as well as the universality of the obtained function 
n g {n) we compare our prediction with experimental data for transverse- momentum 6-jet 
spectra from CMS Collaboration at the CERN LHC Q (Fig. EJ and CDF Collaboration 
at the Fermilab Tevatron [l8j (Fig. [9]). In both cases we find a good agreement between 
theoretical predictions and experimental data. 

Looking in Figs.[5]-|9l we find that contribution of the gluon-to-bottom-pair fragmentation 
in inclusive 6-jet production py-spectra increases from the 10-15 % at the pt — 50 GeV up 
to 30-40 % at the px — 350 GeV. This conclusion contradicts to the prediction of the NLO 
calculations in the collinear parton model in which the gluon-to-bottom-pair fragmentation 
mechanism would be dominant at the large-pr region at the LHC Collider and it would be 



about 50 % at the Tevatron Collider 



26 



4l|. 



Comparing as a whole our results with theoretical predictions obtained in the NLO of 
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a parton model, which also describe ATLAS data for 6-jet production |l4j, we would like 
to pay attention to difficulties of the fixed order collinear calculations. At first, the K- 
factor between LO and NLO calculations is very large at the high p^. The scale uncertainty 
decreases from LO calculation to NLO calculation, but it still remains large. The last one 
can be a signal on large NNLO contributions, which are not taken into account. Second, 
to describe data at non-large px < 50 GeV at the energy of y/S = 2 - 7 TeV in the 
collinear parton model it is needed to add the soft gluon resummation procedure, which is 
far from the application field of DGLAP evolution equation, and should be considered as a 
phenomenological trick rather than rigorous approach. The both these difficulties are solved 
in the PRA by introducing the off-shell LO Reggeized parton amplitudes and unintegrated 
gluon PDF, which take into account large logarithmic contributions in all orders in a s : 
(aM^/^QCD)) n and (a s ln(l/x)) n . 

In such a way, we have obtained the self-coordinated description in the PRA of 66-dijet 
cross sections, where the open 66-quark pair production works solely, and the 6-jet inclusive 
cross sections, where open 66-quark pair production is the main contribution, while the 
gluon-to-bottom-pair fragmentation production is also important. 



IV. CONCLUSIONS 

The CERN LHC is currently probing particle physics at the terascale cm. energies y/S, 
so that the hierarchy A QCD <C \i <C \/S, which defines the MRK and QMRK regimes, is 
satisfied for processes of heavy quark (c or b) production in the central region of rapidity, 
where /x is of order of their transverse momentum. In this paper, we studied QCD processes 
of particular interest, namely inclusive 6-jet and 66-dijet hadroproduction, at LOs in the 
parton Reggeization approach, in which they are mediated by 2 — > 1 and 2 — > 2 partonic 
subprocesses initiated by Reggeized gluon collisions. 

We describe well recent LHC data measured by the ATLAS Collaboration [ljl at the 
whole presented range of the bb-jet transverse momenta, the bb-jet rapidity, the 66-dijet 
invariant mass Mjj, the azimuthal angle between the two jets A0 and the angular variable 



X- We show that the gluon-to-bottom-pair fragmentation component 24j, which takes into 
account effects of large logarithms \og{px / m^) , increases in inclusive b— jet production at 
the high transverse momenta pt up to 30-40 % of sum of all contributions. The extracted 
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by the fit of data the 66-pair multiplicity is in agreement with the previous measurements 



at the LEP Collider 



38 



Comparing different unintegrated gluon PDFs, we have found 



that the agreement with the data has been obtained when we used KMR PDF [31[, and the 



calculations with Blumlein PDF 



35] regularly underestimate data, approximately by factor 



2, in the region of large b— jet pt and the large bb— dijet invariant mass. 
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FIG. 1: The 66-dijet cross-section as a function of dijet invariant mass Mjj for 6-jets with > 



141 ] , the solid polyline corresponds 



40 GeV, \y\ < 2.1. The data are from ATLAS Collaboration 
to KMR unintegrated PDF, the dashed one — to Blumlein PDF. The shaded bands indicate the 
theoretical uncertainties in the case of KMR unintegrated PDF. 
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FIG. 2: The 66-dijet cross-section as a function of the azimuthal angle difference between the two 
jets for 6-jets with pj- > 40 GeV, \y\ < 2.1 and a dijet invariant mass of Mjj < 110 GeV. The data 
are from ATLAS Collaboration 14 ]. the solid polyline corresponds to KMR unintegrated PDF, the 
dashed one — to Bliimlein PDF. The shaded bands indicate the theoretical uncertainties in the 
case of KMR unintegrated PDF. 
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FIG. 3: The 66-dijet cross-section as a function of x f° r ft-jets with px > 40 GeV, \y\ < 2.1 and 



I V boost I 



\ 1 2/1 +2/2 1 < LI, for dijet invariant mass range 110 < Mjj < 370 GeV. The data are from 



ATLAS Collaboration [l4 |. the solid polyline corresponds to KMR unintegrated PDF, the dashed 
one — to Bliimlein PDF. The shaded bands indicate the theoretical uncertainties in the case of 
KMR unintegrated PDF. 
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FIG. 4: The 66-dijet cross-section as a function of x f° r 6-jets with px > 40 GeV, \y\ < 2.1 and 
\Uboost\ = \ I Vi +V2 1 < 1.1, for dijet invariant mass range 370 < Mjj < 850 GeV. The data are from 



ATLAS Collaboration [14]], the solid polyline corresponds to KMR unintegrated PDF, the dashed 
one — to Bliimlein PDF. The shaded bands indicate the theoretical uncertainties in the case of 
KMR unintegrated PDF. 
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FIG. 5: Inclusive differential fo-jet cross-section as a function of px for 6-jets with \y\ < 2.1. The 
data are from ATLAS Collaboration [14]. The dashed polyline corresponds to contribution of the 
open 6-quark production, the dashed-dotted one — the gluon-to-bottom-pair fragmentation, the 
solid — sum of their all. The calculation is done with the KMR unintegrated PDF. 
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FIG. 6: The bb-pair multiplicity n g in a gluon jet as a function of px extracted from the ATLAS data 
for the inclusive 6-jet production spectra The open circles and dashed fitting line correspond to 
Blumlein unintegrated PDF, the black circles and solid fitting line correspond to KMR unintegrated 
PDF. 



21 



> 

CD 

a 

Q. 

X 



CL 

a. 

H 

Q. 
T3 

b~ 

eg 




100 150 200 250 300 350 400 
p T , GeV 



FIG. 7: Inclusive double-differential 6-jet cross-sections as a functions of px for the different rapidity 



ranges: (1) \y\ < 0.3 ( x 10 6 ), (2) 0.3 < \y\ < 0.8 ( x 10 4 ) 
1.2 < \y\ < 2.1. The data are from ATLAS Collaboration 
sum of all contributions (|15|) and KMR unintegrated PDF. 



(3) 0.8 < \y\ < 1.2 ( x 10 2 ) and (4) 



14|. The solid polylines correspond to 
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FIG. 8: Inclusive differential b-jet cross-section as a function of px for 6-jets with \y\ < 2.4. The 
data are from CMS Collaboration 



401 ] . The dashed polyline corresponds to contribution of the 
open 6-quark production, the dashed-dotted one — the gluon-to-bottom-pair fragmentation, the 
solid — sum of their all. The calculation is done with the KMR unintegrated PDF. 
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FIG. 9: Inclusive differential fr-jet cross-section as a function of px for 6-jets with \y\ < 0.7. The 
data are from CDF Collaboration 



181 ] . The dashed polyline corresponds to contribution of the 
open 6-quark production, the dashed-dotted one — the gluon-to-bottom-pair fragmentation, the 
solid — sum of their all. The calculation is done with the KMR unintegrated PDF. 
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